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ABSTRACT 

Hydrodynamic models of short-period Cepheids were computed to determine the pul¬ 
sation period as a function of evolutionary time during the first and third crossings of 
the instability strip. The equations of radiation hydrodynamics and turbulent convec¬ 
tion for radial stellar pulsations were solved with the initial conditions obtained from 
the evolutionary models of population I stars (X = 0.7, Z = 0.02) with masses from 
5.2M 0 to 6.5 Mq and the convective core overshooting parameter 0.1 < a ov < 0.3. In 
Cepheids with period of 4 d the rate of pulsation period change during the first cross¬ 
ing of the instability strip is over fifty times larger than that during the third crossing. 
Polaris is shown to cross the instability strip for the first time and to be the funda¬ 
mental mode pulsator. The best agreement between the predicted and observed rates 
of period change was obtained for the model with mass of 5.4M 0 and the overshooting 
parameter a ov = 0.25. The bolometric luminosity and radius are L = 1.26 • 1O 3 L 0 and 
R = 37.5Rq, respectively. In the HR diagram Polaris is located at the red edge of the 
instability strip. 

Key words: stars: evolution - stars: fundamental parameters - stars: individual: 
a UMi (Polaris) - stars: variables: Cepheids 


1 INTRODUCTION 

Polaris ( a UMi) is a short-period (II sw 3.969 d) small 
amplitude Ce p heid with stead i ly inc r easing pulsation perio d 
llFerniel 1 19661 : lArellano Ferrol Il983l : IDinshaw et al.l Il989h . 
Th e rate of period change is estimated from II = 4 .5 s yr - 1 
dTurner et al.ll2005h to fl = 4.9 s yr -1 (ISprecklev fc Stevensi 
l2008h . These values are over an order of magnitude larger 
than the rates of period change i n most of t h e Ce pheids 
with a period of II ss 4 d. This lead iBerdnikov et al.l rtl997l ) 
to the conclusion that Polaris is undergoing the first cross- 
ing of the ins t ability strip after th e main -sequence phase. 
iNeilson et al.l ll2012f) and iNeilsonl (120141 ') however casted 
doubts on this assumption and argued that Polaris is evolv¬ 
ing along the blue loop. They also showed that the observed 
rate of period change is more consistent with predictions 
from stellar evolution models provided that Polaris is cur¬ 
rently loosing mass at a rate of M ~ 10~ 6 Mq yr -1 . 

Since about the middle twentieth century Po¬ 
laris showed de crease in photometr ic and radial veloc ¬ 
ity amplitudes (lArellano Ferrol Il983l : iFernie et al.l Il993ll . 
IDinshaw et al.l d 1989T) supposed that Polaris is approaching 
the red edge of the instability strip but in 1994-1997 the de - 
crease in amplitude has stopped (Ka mper fe Ferniei ri998f . 
Recent photometric and radial velocity measurements reveal 
slow increase of the pulsation amplitude dBruntt et al.|[2008l : 
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iLee et all 120081 : ISprecklev fc Stevensi 1 20081). The c ause o f 
the amplitude grow t h remains uknown. Fernie et al.l (Il993t) , 
lEvans et al.l (l2002f) . I Turner et al.l (12005 ) reported that Po¬ 
laris lies well inside the instability strip. 

Polaris is the nearest Cepheid and determination of 
its pulsation mode is of great importa nce i n establis hing 
the C epheid period-luminosity relation. iFeast fe Catchpolg 
dl997ll showed that the Hipparcos trigonometrical parallax 
data are consistent with first overtone pulsation in Polaris. 
This conclusion was la ter conf i rmed by nonlinear stellar pul¬ 
sation models by iBono et al.l (l200lh . However recent mea¬ 
surements of photometric and spectroscopic parallaxes seem 
to be more consist ent with fundamental mode pulsation 
(j Turner et al.|[20l3 l. 

The Hip parco s parallax of 7.54 ± 0.11 mas 
(Ivan Leeuwenl l2007l l implies a distance to Polaris of d = 
133 ±2 pc but the spectroscopic and photo metric para l¬ 
l axes give the distance from d = 99 ± 2 pc ( Turner et all 
l2013ll to d = 110 pc llUsenko fe Klochkovall2008l l. There- 
fore the mean angular diameter of 3.28 ma s measured w ith 
the Navy Prototype Optical Interferometer (iNorderen et al.1 
l2000h corresponds to the linear radius of R = 47 Rq for the 
distance d = 133 pc and R = 35Rq for d = 99 pc. 

Polaris is the primary component of a close b inar y sys - 
tem with an orbital period of 29.6 yr. lEvans et al.l (1200811 
directly detected the secondary component using UV im¬ 
ages (A ~ 2255A) obtained with the Hubble Space Telescope 
and found that Polaris has a mass in the range M = 3.1 
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6.7 Mg. This is the first direct estimate of the Cepheid dy¬ 
namical mass. 

In this work, we aim to clarify some contradictions men¬ 
tioned above. We calculate the rates of period change in 
Cepheids with periods near 4 d and determine the evolu¬ 
tionary status of Polaris by comparing the observed rate 
with theoretical predictions. To this end we compute the 
non-linear stellar pulsation models as solution of the equa¬ 
tions of radiation hydrodynamics and turbulent convection 
with initial conditions obtained from the evolutionary com¬ 
putations. Such an approach the author used earlier to ob¬ 
tain theoretical rates of pulsation period change in LMC 
and Galactic Cepheids (lFadevevl2013l . 1201414 . In comparison 
with previous works both evolutionary and hydrodynamical 
computations of this study were carried out with improved 
input physics. As shown below, determination of the evolu¬ 
tionary status provides constraints on the pulsation mode 
and the fundamental parameters of Polaris. 


2 METHODS OF COMPUTATION 
2.1 Stellar evolution 


This work is aimed at determining the pulsation period II 
as a function of evolutionary time t ev during the first and 
third crossings of the Cepheid instability strip. We com¬ 
puted the evolutionary tracks for non-rotating population I 
stars with initial masses 5.2 Mq ^ Mzams ^ 6.5Mg evolv¬ 
ing from the zero age main sequence (ZAMS) until the end 
of core helium burning. Initial mass fractional abundances 
of hydrogen and heavy elements are X = 0.7 and Z = 0.02, 
respectively. The initial relative abundances of the elements 
heavier than heliu m nor malize d by Z were determined ac- 
coring to lAsplund et all l2009lb Evolutionary computations 
were done with a co de implementing the Henyey method 
llFadevevI 120071 . l20ld ) . Below we comment some improve¬ 
ments used in the present work. 

Thermodynamic quantities are obtained by interpola¬ 
tion in th e ta ble s cal culated with the library of programs 
FreeEOS (|lrwinl l2012h . Rosseland opacitie s are computed 
from OPAL data ( Iglesias fe Roeersl Il996l f complemented 
at lo w temperatures by the opacities from (IFerguson et al.l 
l2005fl . Th e nuclear reaction rates are from the NACRE 
database dAngulo et alJ 1 19991 '! with updates for reactions 
14 N(p, 7 ) 15 0, 4 He(aa, 7 ) 12 C and 1 2 C(q, 7 ) 16 Q tha t were 
taken from the REACLIB database dRauscher et al.ll20ld 'l. 

In the evolutionary computations convection is treated 
according to the mixing length theory iBohm-Vitensell 19581 1 
using a mixing-length parameter a a = 1-6. Boundaries 
of convective zones are determined with the Schwarzschild 
criterion but the convective core is extended in radius by 
OotEp, where Hp is the pressure scale height and q ov is the 
dimensionless overshooting parameter ranging from 0.1 to 
0.3. 

To assure the gradual growth of the convective core dur¬ 
ing the core helium burning phase the evolutionary mod¬ 
els were computed with ~ 10 4 mass zones. The mass loss 
ra te M during the main sequen ce phase is from formulae 
b dVink et alJ ( 19991 . l200d . 12001 ) and outside their domain 
M is evaluated according to Nieuwenhuiizen fc de Jagei] 
dl99d) . It should be noted that in our problem the role 


of mass loss is insignificant. For example, the mass of the 
Cepheid Mzams = 5.9 Mq during the first and third cross¬ 
ings of the instability strip is M — 5.89Mg and M = 
5.84Mq, respectively. 

Bearing in mind the importance of initial conditions 
for correct determination of the predicted rates of pulsa¬ 
tion period change we carried out additional evolutionary 
computations for stars with initial composition A' = 0.72, 
Z = 0.014 and the overshooting parameter a ov = 0.1. Re¬ 
sults of these computations were compared with evolution¬ 
ary model s of non-rotating sta rs with initial masses 5 Mq 
and 7 Mq (lEkstrom et alj|2012l4 . The difference in luminos¬ 
ity and time scales was found to be less than 5 per cent. 


2.2 Non-linear stellar pulsation 


To compute non-linear stellar pulsation models we solve the 
Cauchy problem for the equations of one-dimensional radia¬ 
tion hydrodynamics with time-depende nt turbulent convec¬ 
tion treated according to iKuhfui (Il986l ). The equations are 
written as 


d U _ GM r 

dt r 2 

d E dU 

-b P -= 

dt dt 

d Et dV 

~df + Pt ~dF 


1 d 

(p + p.) +[ /„, 
~{r’(F, + F<)]-C„ 


1 d 
pr 2 dr 


(r 2 F t ) +C t + E tv . 


( 1 ) 

( 2 ) 

( 3 ) 


Here t is time, G is the Newtonian constant of gravitation, 
M r is the mass interior to radius r, 17 is the velocity of the 
gas flow, P is the total (gas plus radiation) pressure, p is 
gas density and V = 1/p is specific volume, E is the specific 
internal energy of gas, F r is the radiation flux computed in 
the diffusion approximation. The turbulent pressure Pt is 
related to the specific turbulent kinetic energy Et by Pt = 
\pEt. 

The convective flux is 


Fc = pT n, (4) 

where T is the gas temperature, 

r) Q 

n = ~a A a E H P El /2 ^ = a A a s e\ /2 Cp (V - V a ) (5) 

is the correlation between entropy and velocity fluctuations, 
V = dlnT/dlnP, V a = (dlnT/dlnP) s , Cp is the specific 
heat at constant pressure, S is the specific entropy. 

The turbulent kinetic energy flux is 

Ft = -a A atpHpE 1 t /2 ^. ( 6 ) 

The coupling term is 

Ct = St - Dt - D T , (7) 


where 

St = -^v a n (8) 

tip 

is the source of specific turbulent energy due to buoyancy 
forces, 


Dt 


„ TP 3 / 2 

Qfd -C/^ 

a A Hp 


( 9 ) 
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is the dissipation of specific turbulent energy due to molec¬ 
ular viscosity. The radiative dissipation is 



Trad 


( 10 ) 


where 


ot\ CpkP 2 
~ 7rad 


( 11 ) 


a is the Stefan-Boltzmann constant, k is opacity, g = 
GM r /r 2 is the gravitational acceleration. 

Viscous momentum and energy transfer rates are 


1 d 

4 3 

tdU UY 

r 3 p dr 

^Ptvr 

^ dr r ). 


and 


Etv 


4-Utv(dU_U\ 2 
3 p \ dr r ) 


respectively, where 
Ptv = OLKOLvnpHpE^J 2 


( 12 ) 


(13) 


(14) 


is the kinetic turbulent viscosity. 

Formulae ©, ©, ©, (HU and HU contain dimension¬ 
less parameters oa, a s , at, ad, 7rad and a m . The mixing 
length to pressure scale height ratio is the same as in the 
evolutionary computations: qa = 1.6. P arameters q s , at , 
and aa are set according to IWuchterl fe Feuchtingeil d 1998ll : 


a s = at = 0.6093a s 


0d = 3V 3' 


Appropriate values 7 ra d = 0.4%/3 and a m = 0.1 we re deter ¬ 
mine d from trial computations of Cepheid models dFadevevI 
120141 4 . Thus, the parameter 7 ra d in the expression for the ra¬ 
diative cooling time of convective e lements HU is five ti mes 
smalle r than originally proposed bv IWuchterl fe Feuchtingeil 
(119981') . It should be noted that nearly the same parame¬ 
ters were found t o be most appropria t e for Cepheid mod ¬ 
elling in works by Kollathet_aji d2002l l , ISzabo et aid d2007h , 
ISmolec fe Moskalikl I 20()gT 2010T ). 

In our study we consider the self-excited radial oscilla¬ 
tions. The initial hydrodynamic model is computed by in¬ 
terpolating the zonal quantities of the evolutionary model 
and interpolation errors play the role of small initial per¬ 
turbations. To diminish the amplitude of initial perturba¬ 
tions we used fine zoning in mass in the outer layers of the 
evolutionary model so that as many as 3 ■ 10 3 outer mass 
zones had the relative radius r > 0.1.R, where R is the ra¬ 
dius of the star. All hydrodynamic models were computed 
for N = 500 mass zones. The mass intervals increase ge¬ 
ometrically inward with the ratio of q ss 1.03. The radius 
ri and luminosity L\ at the innermost zone are assumed 
to be time-independent. The radius of the inner boundary 
satisfies the condition r\ ^ O.I.R. 


3 RATES OF PULSATION PERIOD CHANGE 

Together with numerical integration of equations © © 

with respect to time t we calculate the kinetic energy 

N 

E K(t) = ^J2 AM >~ 1 / 2 ^'(*) 2 . ( 15 ) 

i =2 



Figure 1. The pulsation growth rate g as a function of effective 
temperature T e s during the first and third crossings of the in¬ 
stability strip for the evolutionary sequence Mzams = 5.9 Mg, 
Oov = 0.1. Hydrodynamic models are shown by filled circles. 


where A= Mj — Mj _i and Mj is the Lagrangian co¬ 
ordinate of the j—th mass zone. The kinetic energy reaches 
the maximum value -Eiqmax twice per pulsation period II 
and for the initial interval of integration with exponential 
change of I?K,max we evaluate the pulsation growth rate 
g = ndln£ K ,maxM. The period II is determined by a dis¬ 
crete Fourier transform of E^(t). The total sampling interval 
ranges from 10 3 to 10 4 pulsation cycles. 

Fig- [D shows the pulsation growth rate g as a function of 
effective temperature T e s during the first and third crossings 
of the instability strip for hydrodynamic models of the evo¬ 
lutionary sequence Mzams = 5.9M 0 , a ov = 0.1. The evo¬ 
lutionary time t ev , stellar luminosity L and stellar radius R 
at the edges of the instability strip are determined for g = 0 
by linear interpolation between two adjacent hydrodynamic 
models with growth rates of opposite signs. The time spent 
in the instability strip is Af ev ,i = 1.60-10 4 yr during the first 
crossing and At e v ,3 = 4.26- 10 5 yr during the third crossing. 
The ratio At ev 3 /At ev i increases with mass and ranges from 
25 to 30 for 5.6Mq < Mzams < 6.1Mg. 

Evolutionary tracks for stars with initial masses 5.6 Mq 
and 5.9 Mq (a ov = 0.1) are displayed in Fig. [5] The sec¬ 
tions of the evolutionary tracks corresponding to pulsations 
with increasing period (the first and third crossings of the 
instability strip) are shown by dotted lines. Cepheids with 
decreasing periods are beyond the scope of the present study 
and therefore the second crossing of the instability strip is 
not indicated in the plots. 

The upper panel of Fig. [3] shows the temporal de¬ 
pendence of the pulsation period II during the first cross¬ 
ing of the instability strip for the evolutionary sequence 
Mzams = 5.9 Mq, a ov = 0.1. For the sake of convenience 
we use r = t ev — t ev , b as an independent variable, where 
fev.b is the evolutionary time at the blue edge of the insta¬ 
bility strip. Periods of pulsationally stable and pulsationally 
unstable models are shown by open circles and filled circles, 
respectively. 

During the first crossing of the instability strip the star 
begins to pulsate in the first overtone but as the star ap- 
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Figure 2. Evolutionary tracks of stars Mzams — 5.6 Mq and 
Mzams = 5.9 Mq ( (>ov = 0.1). The sections of tracks corre¬ 
sponding to the first and third crossings of the instability strip 
are indicated by dotted lines. 


proaches the red edge it becomes the fundamental-mode 
pulsator. For hydrodynamic models near the pulsational 
mode switching we evaluated the periods of the first over¬ 
tone (lh) and fundamental mode (f) using the additional 
discrete Fourier transform of Ek- The sampling interval was 
limited by the amplitude growth stage when the power spec¬ 
trum reveals the presence of secondary modes with decaying 
amplitudes. Periods of such modes are shown in the upper 
panel of Fig. [3] by open triangles. The boundary between 
the first overtone and fundamental mode pulsators is de¬ 
fined as a middle point between two adjacent hydrodynamic 
models with different orders of the primary mode. For the 
first crossing of the instability strip the second-order al¬ 
gebraic polynomials approximate the temporal dependence 
II(r) with an rms error less than 0.1%. In the upper panel 
of Fig. [3] the analytical fits are shown by solid lines. 

Time derivatives of analytical fits to II(r) are shown in 
the lower panel of Fig. [3] as a function of pulsation period 
II. For the sake of convenience the rate of pulsation period 
change II is expressed in units of seconds per year. Solid 
squares show estimates of II obtained by numerical differ¬ 
entiation. The small deviation of numerical derivatives from 
the analitical fits to II(t) is due to the limited accuracy of 
evolutionary calculations. 

Fig. [4] displays the temporal dependence of II (upper 
panel) and the plot of II as a function of II (lower panel) 
for models of the evolutionary sequence Mzams = 5.9 Mq, 
q ov =0.1 during the third crossing of the instability strip. 
All hydrodynamic models were found to pulsate in the fun¬ 
damental mode. The temporal dependence II(t) is fitted 
with the relative rms error less than 0 . 1 % by the fourth- 
order algebraic polynomial. 

The results of polynomial fitting for hydrodynamic 
models with pulsation periods 2 d ^ II ^ 8 d are displayed 
in Fig. [5] where the rates of period change II are plotted as a 
function of period II for evolutionary sequences with initial 
masses 5.4Mq ^ Mzams ^ 5.9 Mq and the overshooting pa¬ 
rameter a ow =0.1. Each curve in this diagram describes the 
evolution of values II, II as the star crosses the instability 



T (yr) 



Figure 3. The pulsation period II as a function of evolutionary 
time r (upper panel) and the rate of period change n as a function 
of period II (lower panel) during the first crossing of the instability 
strip for the evolutionary sequence Mzams = 5.9 Mq, q ov = 0.1. 
In the upper panel solid lines show the second—order polynomial 
fits, open circles and filled circles show the periods of pulsation- 
ally stable and pulsationally unstable models, respectively. Open 
triangles show the periods of the secondary mode near the mode 
switching. In the lower panel solid lines show time derivatives of 
analytical fits. Solid squares represent the time derivatives ob¬ 
tained by numerical differentiation. 

strip. The solid and dashed lines represent the fundamental 
mode and first overtone, respectively. The plots located in 
the upper part of Fig. [5] represent the first crossing of the 
instability strip and the plots located in the lower part of the 
diagram correspond to the third crossing. As is clearly seen, 
for Cepheids with a period of II ss 4 d the rate of period 
change during the first crossing of the instability strip ITi is 
over fifty times larger than the rate of period change during 
the third crossing II 3 . It should be noted that large values 
of the ratio II 1 /II 3 are also typical for long-period Cepheids 
dFadevevll2014l ). 

The position of Polaris is marked in Fig.[5]for the period 
II* = 3.969 d and the rate of period change II* = 4.7 s yr _1 
which is the mean valu e of the observationa l estim ates by 
iTurner et alJ (120051 1 and lSnrecklev fc Stevensl (120081 1 . Even a 
cursory glance at Fig.[5]leads us to the following conclusions: 
(1) Polaris crosses the instability strip for the first time; (2) 
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Figure 4. Same as Fig. [3] but for the third crossing of the insta¬ 
bility strip. All hydrodynamic models show pulsations in the fun¬ 
damental mode. The solid lines represent the fourth—order poly¬ 
nomial fitting (upper panel) and its time derivative (lower panel). 



logn (d) 

Figure 5. The rate of period change A as a function of pulsation 
period II for evolutionary sequences with q ov = 0.1 during the 
first and third crossings of the instability strip. The numbers at 
the curves indicate the initial stellar mass Mzams hi solar units. 
Solid and dashed lines correspond to the fundamental mode and 
first overtone, respectively. 



-Mzams/M 0 

Figure 6. The ratio of the predicted to observed rates of period 
change TI/TI* as a function of initial stellar mass Mzams for the 
overshooting parameters a ov = 0.1, 0.2, 0.25 and 0.3. 

Polaris is the fundamental mode pulsator. Below we discuss 
in more detail both these conclusions. 

The diagram in Fig.[5]was obtained for the evolutionary 
sequences computed with the overshooting parameter q ov = 
0.1 and the ratio of the predicted to observed rates of pe¬ 
riod change is If/II* « 2. However the rate of period change 
during the first crossing of the instability strip depends on 
the assumed overshooting during the main-sequence evolu¬ 
tionary phase, so that agreement between theory and ob¬ 
servations can be substantially improved. This is illustrated 
by Fig. [3 where the ratio 11/11* is plotted as a function 
of initial mass Mzams for several values of the overshoot¬ 
ing parameter a ov ■ The upper and lo wer errorbars corre¬ 
spond t o the observationa l esti mates bv lTurner et al.1 (120051) 
and ISprecklev fc Stevens! d2008l) . respectively. As seen from 
Fig. [6] the best agreement between the predicted and ob¬ 
served rates of period change (n/n* ss 1.25) is obtained for 
the overshooting parameter a ov = 0.25 and the initial stellar 
mass 5.4 Mq < Mzams ^ 5.45M 0 . 


4 FUNDAMENTAL PARAMETERS OF 
POLARIS 

To determine the pulsation mode of first-crossing Cepheids 
with period of II = 3.969 d we considered the grid of evolu¬ 
tionary and hydrodynamic models with initial stellar masses 
5.2 Mg < Mzams ^ 6.5 M 0 and the overshooting parame¬ 
ter ranging from 0.1 to 0.3. Typical results obtained for the 
overshooting parameter a ov = 0.2 are shown in Fig. [3 

The filled circles connected by dashed lines in Fig. [3 
show the pulsation period at the blue and red edges of the 
instability strip, whereas the open circles connected by dot¬ 
ted lines show the pulsation period at the mode switch. 
Thus, for a fixed value of Mzams the plots in Fig. [3 show 
the range of periods for pulsations in the first overtone and 
in the fundamental mode. The similar diagrams were ob¬ 
tained for other values of the overshooting parameter a ov 
and no models pulsating in the first overtone with period of 
n = 3.969 d were found. Therefore, Polaris is undoubtedly 
a fundamental-mode pulsator. 
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Figure 7. The pulsation period at the blue edge (lower dashed 
line) and red edge (upper dashed line) of the instability strip for 
models of first—crossing Cepheids computed with the overshoot¬ 
ing parameter q ov = 0.2. Lower and upper dotted lines show the 
periods of the first overtone and fundamental mode at the mode 
switch boundary. Filled and open circles represent the interpo¬ 
lated values of the period. 


Figure 8. Evolutionary tracks of stars with initial masses 5.5 Mq 
and 6.2 M® (a ov = 0.2) near the first crossing of the instability 
strip. Dashed lines show the blue and red edges of the instability 
strip. The dotted line shows the mode switch boundary. The star 
symbols show the location of models with the pulsation period 
II* = 3.969 day 


Table 1. The upper and lower values of the fundamental param¬ 
eters of Polaris 


Qov 

M/M® 

L/L e 

R/R® 

T e a 

(K) 

d 

(pc) 

0.1 

6.38 

1770 

39.9 

5930 

113 


6.08 

1390 

39.3 

5630 

111 

0.2 

5.88 

1660 

38.7 

5930 

110 


5.59 

1300 

37.8 

5640 

107 

0.25 

5.59 

1510 

37.9 

5850 

107 


5.39 

1260 

37.4 

5640 

106 

0.3 

5.39 

1450 

37.5 

5820 

106 


5.19 

1220 

36.7 

5630 

104 


Intersection of the dashed and dotted lines with the 
solid horizontal line II = 3.969 d in Fig. [7] gives the lower 
and upper values for the initial stellar mass Mzams, respec¬ 
tively. The upper and lower values of the mass M, radius 
R and luminosity L are evaluated by interpolation between 
the evolutionary tracks. 

Fig.[8]shows the evolutionary tracks in the HR diagram 
for stars with initial masses 5.5 M® and 6.2 M® (q ov = 0.2). 
The dashed lines show the blue and red edges of the insta¬ 
bility strip. The mode switch boundary shown by the dotted 
line is not the straight line due to difficulties in determin¬ 
ing the exact location of the point of mode switch. Models 
with pulsation period II* = 3.969 day are shown by the star 
symbol. 

The lower and upper values of fundamental parameters 
for Polaris are given in Table [I] The last column gives the 
distance to Polaris d ev aluated from the obse rved angular 
diameter 6l = 3.28 mas (iNordgren et alJfeOQOl ). 


5 CONCLUSION 

The evolutionary and pulsational models together with the 
observed period n* allow us to obtain the lower and upper 
estimates of the initial mass of Polaris. As seen in Fig. [6] a 
major cause of uncertainty in the determination of mass is 
the overshooting parameter a ov which is the unknown a pri¬ 
ori quantity. The comparison of theoretical predictions with 
the observed rate of period increase leads to the narrower 
mass range 5.2 M® Mzams ^ 5.7Mq and to the constraint 
on the overshooting parameter a ov ^ 0.2. The best agree¬ 
ment with observations (n/II* ~ 1.25) was found for the 
model Mzams = 5.4Mq, a ov = 0.25 which locates near the 
red edge of the instability strip. The mass, M = 5.39 M®, 
is within the range 3. 1M® ^ M ^ 6.7Mq of the measured 
dynamical mass (lEvans et al.ll2008l ). 

The mean linear radius R = 37.5 R® for the overshoot¬ 
ing parameter a ov = 0.25 (see Table [lj and t h e observed 
angular diameter 9l = 3.28 mas (iNorderen et al.1 hood ) 
lead to the distance of d = 106 pc. This value is consis¬ 
tent with photometric and spectroscopic parallaxes imply- 
ing a distance from 99 pc (ITurner et al.l 120131) to 110 pc 
(lUsenko fc Klochkova! 120081) . To further improve theoreti¬ 
cal predictions one has to take into account effects of ini¬ 
tial chemical composition. For example, preliminary com¬ 
putations showed that the increase of the fractional mass 
abundance of hydrogen from A' = 0.7 to X = 0.73 leads to 
decrease in IT by « 12%. 

The fact that Polaris is evolving on a Kelvin-Helmholtz 
timescale after the main sequence phase and is crossing the 
instability strip for the first time allows us to conclude that 
effects of mass loss do not play a significant role in evolu¬ 
tionary changes of the pulsation period. 

Decrease of the ratio IT/IT* with decreasing initial mass 
Mzams at the fixed value of the overshooting parameter 
favours the closer location of Polaris to the red edge of 
the instability strip. This assumption does not necessarily 
contradict the growth of the pulsation amplitude observed 
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for two last decades. Indeed, the balance between the k- 
mechanism exciting pulsational instability and convection 
which suppresses oscillations can cyclically change due to 
variations of magnetic field on time scale longer than the 
pulsation period. Therefore, the currently observed increase 
of the pulsation amplitude might be due to the gr owth of a 
magn etic field which tends to weaken convection dStothersI 
120091 b Longitudinal magnetic fields from 10 to 30 Gauss 
were detected in P olaris by llBorra et al.l 1 198 lib whereas 
lUsenko et al.l ([2010) reported on detection of variable mag¬ 
netic field from -109 to 162 Gauss. 

T he present study is based on the convection model 
by iKuhfufj (Il986h which contain several parameters of or¬ 
der unity. The most uncertain of them are 7 ra d and a m in 
expressions for the radiative cooling time CP and kinetic 
turbulent viscosity HP, respectively. The values of these pa¬ 
rameters affect the damping of oscillations and are inferred 
from observational constraints (for example, the width of 
the instability strip). At the same time, the uncertainty in 
the values of 7 ra d and a m does not play an important role in 
determinating the pulsation period. This is due to the fact 
that the depth of the convective zone in Cepheids is only a 
few per cent of the stellar radius, whereas the length of the 
fundamental mode period is of the order of two sound travel 
times between the centre and the surface of the star. 
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